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Abstract 

We compare the behavior of overlap fermions, which are chirally invariant, 
and of Wilson twisted mass fermions at full twist in the approach to the chiral 
limit. Our quenched simulations reveal that with both formulations of lattice 
fermions pion masses of 0(250 MeV) can be reached in practical applications. 
Our comparison is done at a fixed value of the lattice spacing a ~ 0.123 fm. 
A number of quantities are measured such as hadron masses, pseudoscalar 
decay constants and quark masses obtained from Ward identities. We also 
determine the axial vector renormalization constants in the case of overlap 
fermions. 
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1 Introduction 



Approaching the so-called physical point in lattice QCD, where the light meson 
masses assume their values observed in nature, is a major challenge for numerical 
simulations in lattice field theory. In this regime of QCD, where the light quarks, i.e. 
the u and d flavors, play the dominant role, the explicit breaking of chirality on the 
lattice, the role of topological charge excitations and the appearance of unphysical 
very-low lying eigenvalues of the lattice Dirac operator can render simulations very 
demanding, if not impossible. 

The complicated low energy dynamics of light quarks of the fundamental theory 
can be effectively described by chiral perturbation theory jTj. In chiral perturbation 
theory (x?T) an effective Lagrangian is constructed with terms that are compatible 
with the global symmetries of the original QCD Lagrangian and that are ordered 
according to some hierarchy, which depends on the expansion regime. In this effec- 
tive Lagrangian, each term has a coefficient which enters as a free parameter that 
cannot be determined within the chiral perturbation theory itself. Such coefficients 
are denoted as the low energy constants, and they play a vital role in many QCD 
processes at low energy. The challenge is then to determine these low energy con- 
stants from first principles, i.e. directly from QCD as the underlying, fundamental 
theory. 

Such a link between the effective theory and fundamental QCD can, in principle, 
be provided by lattice techniques. However, simulations at values of quark masses, 
where such a contact to ;\PT can safely be made, are very difficult, as mentioned 
above. Standard approaches using the Wilson fermion discretization of lattice QCD 
are confronted with technical algorithmic problems since in this case the bare quark 
mass does not provide an infrared regulator. In this paper we will consider two 
formulations of lattice QCD that are able to overcome this problem: overlap fermions 
|2] and twisted mass fermions [3]. In particular, we will concentrate on the Wilson 
twisted mass (Wtm) [3] formulation of lattice fermions at a full twisting angle of tt/2. 
Both, overlap and Wtm fermions lead to 0(a) improvement for physical quantities 
and are expected to allow for simulations at very small quark masses, corresponding 
to their physical values as estimated from experiment. Overlap fermions induce an 
exact lattice chiral symmetry j^j and provide a sound definition of the topological 
charge Hence they have conceptual advantages, but they are, unfortunately, 
rather expensive in numerical simulations. Twisted mass fermions on the other hand 
are rather cheap to simulate but show residual chiral symmetry breaking effects. 
Hence it is not clear how the benefits of both formulations of lattice QCD compare 
in practical simulations and here we would like to give a first direct comparison. 

We compute a number of quantities, Ward identity quark masses, meson and 
baryon masses, decay constants and renormalization factors, driving the values of 
quark masses small enough that xPT is expected to be applicable. While for Wtm 
fermions this is the first work along this line, we refer to Refs. [7j for other simulations 
in the p-regime of chiral perturbation theory in the case of overlap fermions. 
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We emphasize that the present comparison between both lattice fermions is only 
performed in the quenched approximation at (3 = 5.85, corresponding to a lattice 
spacing of a ~ 0.123 fm. No attempt of a scaling analysis is performed here, see 
however Ref. [H] for a first work in this direction for Wtm. It is the main goal of this 
paper to investigate how both formulations of lattice QCD behave in their approach 
to the chiral limit. In particular we are aiming at an investigation to how small 
quark masses both formulations can be driven when used in practical simulations. 

We will also provide a computing time estimate from our results in Ref. jHj- This 
question is most important for eventual dynamical simulations. If, for many quan- 
tities, Wilson twisted mass fermions can reach quark masses that can be compared 
with overlap fermions, their advantage in the simulation cost will immediately ben- 
efit in dynamical simulations. For first results of dynamical twisted mass fermions, 
see Ref. [TU]. 

2 Overlap and Wilson twisted mass fermions 

In this section we discuss the two fermion formulations of lattice QCD that we have 
employed to study the approach to the chiral limit at small values of the quark mass: 
overlap fermions and Wilson twisted mass fermions. We emphasize some relevant 
points for the present work and refer to Refs. Ej for further discussions. 

2.1 Overlap fermions 

Over the last few years a lattice fermion formulation leading to an exact lattice 
chiral symmetry was elaborated, namely the Ginsparg- Wilson fermions. The cor- 
responding lattice Dirac operator Dqw satisfies the Ginsparg- Wilson relation ^3] 

£>gw75 + lbD G w = 2aD GW >y 5 RD GW , (1) 

where R is a local term. The realization of an operator -Dgw that we use here is the 
overlap fermion, which is characterized by the Neuberger-Dirac operator |2J. For 
Rxy — $xyl (2p) it takes the form 

D ov = (l-I^)D(°)+™ ov , 

D<$ = ^{l + A/ViU}, A = aD w -p, (2) 

where a = a/p and Dyj is the standard Wilson-Dirac operator, 

4 i 

D * = Y, + v M ) - av;v,]. (3) 

V M and V* denote the usual forward and backward covariant lattice derivatives, 
m ov is the bare quark mass and p > 1 is a mass parameter, which we set to 1.6. 
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Dov (the overlap operator at zero quark mass) obeys the Ginsparg- Wilson rela- 
tion in Eq. (JTJ. It does obey a lattice modified but exact chiral symmetry, which 
turns into the standard chiral symmetry in the continuum limit [5]. This symmetry 
protects the lattice fermion from additive mass renormalization and from 0(a) lat- 
tice artifacts, i.e. the action built by the operator in Eq. (J2J) is 0(a) improved. It 
also implies that there are exact zero modes with a definite chirality Thus the 
topological charge can be identified as the index obtained from these zero modes. 
0(a) improved bilinears are constructed as follows: 

o-fr(i-^y = (rr^) (4) 

where ip a and ip 13 are two different flavours, and the last equality holds for correlation 
functions at non-zero physical distance. In the following we will use the following 
notation: 0°; = P, Of = S, O™ = V„ and O^ = A M . 

By now the overlap fermion has a very well established theoretical basis, but its 
simulation is rather tedious. In our code the inverse square root is approximated by 
Chebyshev polynomials to an absolute accuracy of 10 -15 (see Ref. jH] for details). On 
this level of precision the lattice chiral symmetry is certainly reliable. Unfortunately 
the computational effort exceeds the one for Wtm by a large amount, which means 
that, at least for large volumes, for the time being only quenched QCD simulations 
are possible with Ginsparg- Wilson fermions, at least when chiral symmetry is to 
be realized to the precision enforced in this paper. The virtues of this formulation 
include also a protection against exceptional configurations, thus allowing simula- 
tions at small pion masses, certainly when they are comparable or larger than their 
physical values. 

2.2 Wilson twisted mass fermions 

As an alternative to regulate exceptionally small eigenvalues, Wilson fermions with a 
twisted mass E] can be used. This means that the Wilson-Dirac operator obtains 
a mass term of the form 

m tm + ifi^r 3 , (5) 

where m tm is again the bare quark mass, r b are the usual Pauli matrices acting in 
flavor space and fi is the "twisted mass" . 

In this paper we will work with Wilson twisted mass fermions that can be ar- 
ranged to be 0{a) improved without additional improvement terms. To be more 
precise, let us start by writing the Wtm QCD action (in the twisted basis) as 

Sp.ijj^] = a A ^2i>(x)(D w + m tm + z/i7 5 r 3 )^(x) , (6) 

X 

where the operator D\y is given in Eq. Q. 
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The action as it stands in Eq. (jSJ can, of course, be studied in the full parameter 
space (m tm , /i). A special case arises, however, when m tm is tuned towards a critical 
bare quark mass m c . In this, and only in this situation, all physical quantities are 
(or can easily be) 0(a) improved. It is hence natural to rewrite 

m tin = m c + m (7) 

with m an offset quark mass. The values of m c need only to be known with 0(a) 
accuracy j3] and can be, for instance, taken from the pure Wilson theory at fi = 0. 

In this approach there is no need of improving the operators, so we will consider 
the usual local bilinears 

P b = ^ 75 y^ S b = ^ 

4 = ^75^ Vl = (8) 

where b is a SU (2) flavour index. 

Of particular interest is the PCVC relation. In the twisted basis it takes the 
form 

dX = -2/ie 36c P c , (9) 

where <9* is the lattice backward derivative. Through a vector variation of the 
action one obtains the point-split vector current as defined in Ref. [3]. This current 
is protected against renormalization and using the point-split vector current, the 
PCVC relation is an exact lattice identity. This implies that Zp = Z' 1 , where Z^ is 
the renormalization constant for the twisted mass fi. This will become important in 
the extraction of the pseudoscalar decay constant f v as described below. Recent tests 
have revealed a very promising scaling behavior of this lattice fermion formulation 
(in quenched simulations) jH]. 



3 Numerical results 

Using standard heat bath and over-relaxation techniques to generate gauge field 
configurations with the Wilson plaquette action, we have performed various simula- 
tions at (3 — 5.85, which corresponds to a value of the lattice spacing a ~ 0.123 fm 
(a -1 ~ 1.605 GeV using r = 0.5 fm [Ej). Periodic boundary conditions were used 
in this work. 

For overlap fermions we have 140 configurations on a 12 3 x 24 lattice (L 12 ~ 
1.48 fm). The bare quark masses are m ov a = 0.01,0.02,0.04,0.06,0.08,0.10. The 
simulations for twisted mass fermions are done at full twist. This is achieved by 
choosing m tm = m c in the pure Wilson theory at fi = 0. In the standard hopping 
parameter k notation the offset quark mass is related to k by fh = ^: — ^r- Choosing 
= then corresponds to use k — k c — 0.16166(2) from the vanishing of 
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the pion mass for Wilson fermions |Hj. The twisted quark mass parameter was 
chosen to be \m = 0.005,0.01,0.02,0.04,0.06,0.08,0.10. We then accumulated 140 
configurations on 12 3 x 24, 140 configurations on 14 3 x 32 (Lu ~ 1.72 fm) and 380 
configurations on 16 3 x 32 (L 16 ~ 1.97 fm) lattices. In tables and plots below both 
m ov and /i are called mb are - 

Our simulations, performed at a number of bare quark mass values, applied a 
multiple mass solver (MMS) in both cases. The MMS for twisted mass fermions will 
be described in a forthcoming paper whereas the MMS for the usual fermion 
formulations can be found in Ref. ^H]. In order to improve the projection on the 
fundamental states of our correlators, we have implemented the Jacobi smearing in 
the way it is described in Refs. [TT)| ITTj. 

In the following we will extract several physical quantities (hadron masses and 
decay constants) using the overlap and the Wtm formulation of lattice QCD. Since 
the flavor structure of the quark propagators is important in the Wtm case, we will 
specify the flavor content of the local operators only for this formulation. We will 
indicate the fermionic action used in order to distinguish the correlation functions, 
where necessary, of the two formulations. 



3.1 Meson masses 

The first quantities we compared are the meson masses M meson . They are extracted 
by fitting the suitable two-point correlation functions to the standard expression 
obtained from a spectral decomposition and retaining only the fundamental state 
(a < x < T) 



a 3 £(Ot(x,* )O(0)> = K0|Q^)l 2 e -M_| CQsh 



-^mes 



T 

-^meson ( %0 



(10) 



in the time interval [t m j n , -|] (with a <C £ min < -|). £ min has been chosen by con- 
sidering the effective mass, the dependence of the fit on t min and by comparing 
with a two-state fit. We extract the pseudoscalar masses (both degenerate and 
non-degenerate) from the correlation functions 

C b PMl (x ) = a 3 J2( pb & x o)P\0)) tm , 6 = 1,2. (11) 

X 

For overlap fermions we consider 

Cp,ov(:co) = a ^(P t (f,x )P(0)) ov , (12) 

X 

and 

Cp-s,ov(xo) = a 3 ^(P f (f , xo)P(0) - S\x, x )S(0)) ov , (13) 



X 
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M p n a 

7T,OV 


M p ~ s a 

7T,OV 


Mif m x24 a 

7r,tm 


7T,tm 


7r,tm 


7T,tm -LD 


0.005 










0.1700(25) 


2.7 


0.01 


0.212(9) 


0.140(20) 


0.2327(70) 


0.2301(37) 


0.2254(19) 


3.6 


0.02 


0.237(7) 


0.196(14) 


0.3193(48) 


0.3175(30) 


0.3122(16) 


5.0 


0.04 


0.299(5) 


0.280(10) 


0.4520(40) 


0.4506(23) 


0.4452(14) 


7.1 


0.06 


0.355(4) 


0.346(8) 


0.5596(35) 


0.5575(19) 


0.5535(12) 


8.9 


0.08 


0.405(4) 


0.401(7) 


0.6541(31) 


0.6510(17) 


0.6488(11) 


10.4 


0.10 


0.450(4) 


0.451(6) 


0.7417(26) 


0.7378(16) 


0.7359(11) 


11.8 



Table 1: We show numerical results for the pion masses obtained with overlap 
fermions using both the correlators Cp jOV and, Cp_s,ov, and with Wtm fermions us- 
ing three different volumes. The overlap results were obtained on a 12 3 x 24 lattice, 
and we set (3 = 5.85 everywhere. The bare quark mass mb are corresponds to m ov in 
the overlap and to \i in the twisted mass case. 

where the contribution of the topological zero modes cancels |18| Ej. This last 
method has the drawback that the scalar meson appears as an excited state and can 
affect the extraction of the mass of the ground state for large quark masses. In our 
study, however, the quark masses are sufficiently small such that this problem never 
occurs. 

In Figure^and Tables^lSlwe give an overview of our results for the pseudoscalar 
meson masses. For twisted mass fermions we have performed simulations on three 
volumes, such that we can investigate finite volume effects. On the 12 3 x 24 lattice, 
for the largest bare quark masses there are small contaminations from the excited 
states (which include also states of opposite parity appearing as 0(a 2 ) artifacts) 
which are removed by using sink-smeared correlators. Apart from this effect, small 
finite volume effects are visible at the smallest masses. An analysis along the lines 
of Ref. ^1 shows that, on the pion masses corresponding to smallest value of M n L 
(i.e. those corresponding to \ia = 0.01 for L = L i2 and to jia = 0.005 for L = L w ), 
the finite volume effects are within 2-3 percent and at most within two standard 



^bare*^ 






M 12 " x24 a 

JW 7r,tm u 


M 14 " x32 a 

Jla 7r,tm u 


M 16 " x32 a 

JW 7r,tm u 


0.01 


0.134(22) 


1.6 


0.2332(70) 


0.2303(39) 


0.2257(21) 


0.02 


0.192(16) 


2.3 


0.3203(49) 


0.3172(28) 


0.3126(18) 


0.04 


0.275(12) 


3.4 


0.4523(41) 


0.4493(22) 


0.4455(16) 


0.06 


0.342(10) 


4.2 


0.5584(38) 


0.5564(20) 


0.5538(14) 


0.08 


0.397(10) 


4.8 


0.6510(36) 


0.6504(19) 


0.6492(13) 


0.10 


0.448(8) 


5.4 


0.7353(35) 


0.7363(18) 


0.7361(13) 



Table 2: The same as Tabled but with sink smearing. 
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deviations from the extrapolated infinite volume limit. In practice they are thus not 
really relevant for the following discussion. Below, we will present only those Wtm 
results which were obtained on the 16 3 x 32 lattice, for which finite volume effects 
are completely negligible as long as \xa is larger than 0.005. 

The problem of isolating the ground state does not appear in the overlap case 
when we analyze the Cp iOV (xo) correlator. Even when we introduce a parity violating 
term explicitly with the C , p_s iOV (xo) correlator, we still do not have any problem to 
extract the pion mass for all the values of the quark masses simulated. This is due 
to the fact that, at a fixed value of the bare quark mass, the pseudoscalar masses 
are smaller (and the gap between the fundamental and the excited states in the 
pseudoscalar correlator or between the fundamental state of the pseudoscalar and 
of the scalar correlators larger) than the corresponding masses for Wtm. 

Given the value of M n L in the case of overlap fermions and the experience from 
Wtm fermions, we expect very small finite volume effects for the 5 heaviest quark 
masses (at the level of few percents in the case of m ov a = 0.02). For the lowest 
quark mass (for which M ntOV L = 1.6) finite volume effects can be more relevant and 
thus we usually do not include the corresponding data point in the fits. However, 
the analysis of the various quantities presented below suggests that, for this quark 
mass, finite volume effects are not larger than our statistical error. 

Note that for 0(a) improved Wilson fermions results j2Hl I2*T] . which are also 
shown in Figure ^ the simulations had to be stopped at rather large values of the 
quark mass to avoid the appearance of exceptional configurations. On the contrary, 
with both Wilson twisted mass fermions and overlap fermions we can reach very low 
values of the quark and hence of the pion mass. 

For overlap fermions, the masses extracted from Cp_s jOV (a;o) have, to a very good 
approximation, a linear behavior down to the smallest mass (M n ~ 220 MeV). A 
linear extrapolation to the chiral limit gives an intercept of —0.002(6). For Wilson 
twisted mass fermions we have performed two fits: a linear one on the four smallest 
masses from which we get an intercept of 0.0017(2) with \ 2 /d.o.f. = 1.35, and a 
quadratic one on all of the 7 points from which we get an intercept of 0.0045(4) 
with x 2 '/d.o.f. = 0.19. As we can see, the data show a behavior which is much 
better described by a parabola than by a straight line. We attribute the value of 
the intercept, non-compatible with zero in the chiral limit, to the 0(a) error in k c , 
which at = gives an 0(a) residual pion mass. The dependence of M 2 upon ji 
and a at small values of the twisted mass can in principle be computed in x?T with 
the inclusion of the twisted mass term [22] • The curvature of the Wtm data at high 
values of /i, absent for overlap fermions, can be explained by the uncertainty of 0(a) 
in k c (determined in the pure Wilson case /i = 0). This induces an uncertainty of 
0(a 2 fi 2 ) in the pion mass, which increases with a higher value of the twisted mass. 

We clearly observe that the pion masses obtained with twisted mass fermions 
shown in Figure Q always lay above the ones from the simulations using overlap 
fermions. This suggests that the renormalization factor Z m of the quark mass should 
be larger for Wilson twisted mass fermions when compared to overlap fermions. It 
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Figure 1: Comparison of quenched results for the pion mass squared as a function 
of the bare quark mass for three lattice fermions: standard 0(a) improved Wilson 
fermions \2IA \21\j . twisted mass fermions and overlap fermions. The bare quark mass 
corresponds to m ov in the overlap and to fi in the twisted mass case. 
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Figure 2: Twisted mass fermions: linear and quadratic fit of the squared pion masses 
as function of the bare quark mass. 
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ivi p a 


M tm n 
ivi p u 


o 005 




U.OuU 1 4:0 I 


0.01 


0.632(34) 


0.468(46) 


0.02 


0.638(26) 


0.543(24) 


0.04 


0.653(16) 


0.6482(91) 


0.06 


0.666(12) 


0.7333(57) 


0.08 


0.683(9) 


0.8087(41) 


0.10 


0.702(8) 


0.8799(33) 



Table 3: Vector meson masses with sink smearing. 

also means that smaller values of the quark mass have to be simulated with Wilson 
twisted mass fermions to reach the same pion mass as with overlap fermions. 

As a next quantity we consider the vector meson mass. It has been extracted 
from the following correlators, 

k=l x 

a 3 3 

CV,ov(x ) = y^T]T(^(f,x )V fc (0)> ov . 

k=l x 

In order to extract a reliable value for the vector meson mass, sink-smearing has 
been used. For overlap data, where the volume is not so large and the statistical 
errors are consequently significant, the smearing procedure is important in order 
to isolate the ground state (at the same time reducing the statistical error of the 
estimated mass). For twisted mass data on the largest volume (16 3 x 32), smearing 
has a small effect, apparently because it is not able to decrease substantially the 
coupling of the interpolating operators to the excited states. The vector meson 
mass has been extracted by averaging the values of the effective mass within a time 
interval which excludes the points at the largest times. This points are in fact 
affected by large statistical uncertainty and tend to lower the estimated value. The 
results are reported in Table El and plotted in Figure El 

In the case of Wtm fermions, we report in Fig.0]two examples for the plateaux of 
the effective mass of the vector meson, namely the case corresponding to the smallest 
quark mass \xa = 0.005 (which has highest statistical fluctuations and for which we 
have chosen the plateau in the range 8 < xq < 13) and that of an intermediate mass 
\ia = 0.06 (with plateau chosen in the range 9 < xq < 12). Despite the precautions 
described above, we still observe a strong decrease of the vector meson mass at 
low quark masses, as compared to the overlap case where the behavior with the 
pion mass squared is linear to a good approximation. We will discuss later, when 
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Figure 3: Comparison of results for the vector meson mass as a function of the pion 
mass squared for overlap and twisted mass fermions. The chiral extrapolation in the 
overlap case is done with a linear fit in M%. 

examining f n , the problems that can arise in the Wtm case when going down to low 
masses at fixed lattice spacing. This phenomenon has to be further investigated, 
for example by going to higher f3 values. This is beyond the scope of the present 
work and will be addressed in another (presently on-going) project whose ultimate 
goal is the study of the scaling behavior of various quantities in twisted mass QCD 
computed for a wide range of quark masses |23j . 

3.2 Renormalization constants 

Using the method explained in Refs. (2011^ we computed the renormalization group 
invariant (RGI) quark mass renormalization constant Z^ 1 for Wtm and overlap 
fermions. Here we just describe the method and we refer to Refs. |2U[ |2"T] for a 
more detailed explanation. There are essentially two requirements for appling this 
method: 

• Due to the symmetries of the lattice action, renormalization constants of dif- 
ferent local operators are related in a simple way. 

• Using an alternative discretization of QCD it is possible to compute a renor- 
malized matrix element (or quark mass) in the continuum. 

The basic idea is to compute a universal factor — that could be a RGI matrix el- 
ement or quark mass — in the continuum at a fixed reference value of a physical 
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Figure 4: Effective mass plateaux for the vector meson mass at jja = 0.005 and at 
fia = 0.06. Notice the different scale on the y axis. 



quantity (for example the pion mass). Then matching this universal factor with a 
renormalized matrix element (the bare matrix element computed in the target regu- 
larization) at the reference point and at fixed lattice spacing, it is trivial to extract 
the renormalization factor. Obviously the method works, using the symmetry prop- 
erties of the target regularization, if it is possible to relate renormalization factors 
of different local operators. To this end, we have used two different matching con- 
ditions. We have matched the RGI quark mass (method 1) and the matrix element 
of the pseudoscalar density (method 2) (see j2I] for details) at the reference points 
given by x rc f = (r M n ) 2 = 1.5736, 3.0, 5.0, where the last point is considered only 
for twisted mass fermions (in the overlap case we do not have data in the region cor- 
responding to such high pion masses). The universal factor was obtained using the 
renormalization constants and the bare matrix elements computed by the ALPHA 
collaboration j2UES3 using 0(a) non-perturbatively improved Wilson fermions. For 
overlap fermions, due to the lattice chiral symmetry j^J, we have the following rela- 
tion between the renormalization factors of the pseudoscalar and scalar density and 
the renormalization factor of the quark mass, 

Z P = Z s = -L . (14) 

For Wtm, due to the existence of an exact flavor symmetry for massless Wilson 
quarks, we obtain 

Zp = ±- (15) 

for all the flavor components. We summarize our results in Table |3] For overlap 
fermions, these results are basically independent from choosing x rc { = 1.5736, 3.0. In 
the twisted mass case, the slight dependence (always within 1-2 standard deviations) 
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^ref 


method 1 

^RGI,ov 
m 


method 2 
7-rgi.ov 


method 1 

yRGI.tm 


method 2 

yRGLtm 


1.5736 


1.02(6) 


0.98(5) 


2.27(7) 


2.22(8) 


3.0 


0.98(7) 


1.01(5) 


2.32(6) 


2.36(6) 


5.0 






2.39(5) 


2.55(20) 



Table 4: Results for Z* GI ' ov and Z* GI ' 





m PCAC a 




0.005 




0.008303(7) 


0.01 


0.00695(2) 


0.016602(8) 


0.02 


0.01391(3) 


0.033187(17) 


0.04 


0.02795(5) 


0.066514(24) 


0.06 


0.04218(6) 


0.100051(35) 


0.08 


0.05659(7) 


0.133899(42) 


0.10 


0.07116(8) 


0.168110(51) 



Table 5: Ward identities quark masses. 



is due to the C(a 2 /x 2 ) lattice artifacts that affect the quantities used for the matching. 
The large error on ,tm at x ro f = 5.0 from method 2 comes from the error on 
the universal factor computed from the data of the ALPHA collaboration in the 
continuum limit. The rather large value of the renormalization constant in the case 
of Wtm fermions is reflected in the slope of the curves shown in Figures ^ and 121 



3.3 Ward identities quark masses 

The Ward identities (WI) quark masses rrip CAC anc ^ 

be extracted from 

the ratios 



m PCAC 



E £ (doAl(x,x ) P(0)> 



2E,(Pt(x> )P(0)) ' 

^ c E.W(x,x ) no)) 

mpcvc " 2£,(p C (f,* )p C (0)> • (16) 

Results for the (WI) quark masses are reported in Table El and plotted in Figure El 
together with the results of a quadratic extrapolation to the chiral limit. 

From the values of the WI quark masses and the bare quark masses, the corre- 
sponding renormalization factors can be computed at each value of the bare quark 
mass. They are then given by 

gov _ m ov 2 tm — (IJ) 

'"PCAC '"PCVC 
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Figure 5: Comparison of results for the WI quark masses as a function of the quark 
mass for overlap and twisted mass fermions. The values of m PC Ac an d rapcvc 
extrapolated quadratically to the chiral limit are 0.00004(2) and 0.00003(1) for the 
overlap and the Wtm case respectively. 

The behavior of the WI quark masses in Figure 03 is only apparently linear. In the 
overlap case this behavior is well described by a quadratic curve (where the quadratic 
term is quite small). This is reflected in the linear behavior of Z°^ with a rather 
mild slope (see Figure E}. A linear fit (excluding the data point at am ov = 0.01) 
allows then to obtain a value for Z°^ in the chiral limit, Z°^ = 1.448(4). In the 
twisted mass case even a quadratic fit is inadequate to describe the behavior of the 
PCVC mass with respect to the bare mass \i. This is reflected in the behavior of 
Zy™ as function of /i — see Figured — which, in particular at small quark masses, 
is altered by residual lattice artifacts. This affects the denominator of Eq. (|17|) (i.e. 
m pcvc)' gi vm g rise to a strong non-linearity. For this reason we prefer not to quote 
any value for Zy^ in the chiral limit. 

3.4 Pseudoscalar decay constants 

The pseudoscalar decay constants can be computed from the ratios 

= g^jojA^ry tm = ^ri(oin & k)t m i b=12 (18) 

7T 7T 

which require the determination of Z™ and of Zy 1 as discussed in the previous 
section (where at maximal twist Zy 1 should be identical to the value computed 
with standard Wilson fermions). 
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Figure 6: as function of the quark mass and its chiral extrapolation. 
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Figure 7: Z^ 1 as function of the quark mass. 
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Figure 8: Comparison of results for the pion decay constant as a function of the pion 
mass squared for overlap and twisted mass fermions. The vertical line represents the 
value of (M n a) 2 which corresponds to a bare quark mass rh q a = a 2 Aq CD; i.e. the 
r.h.s. of Eq. \21\) apart from an unknown proportionality factor. 

There is also a second, "indirect" method that uses the PCAC and PCVC rela- 
tions and does not require the computation of any renormalization constant, 

^=^K0|Pk) ov |, r = ^l(0|PV)t m |, 6 = 1,2. (19) 

This second method will prove to be useful especially for the Wtm case. 

From the theoretical side we expect from one loop quenched chiral perturbation 
theory (q%PT) that f n has neither chiral logarithms nor finite volume effects, 

u=s i} + jSjf M ^) • < 20) 

where a 5 is a low energy constant. 

As we have seen in the previous section, a reliable number for Z^ 1 cannot be 
provided given the strong quark mass dependence of the renormalization constant 
in the Wtm case. We therefore only give values for /„. from the "indirect" method. 
Results for the pion decay constant obtained with the "indirect" method for both 
overlap (from Cp_s,ov(xo)) and twisted mass fermions (from Cp tm (x )) are reported 
in Table lOand plotted in Figure El 

This Figure reveals that in the case of the overlap fermion f n nicely follows the 
linear behavior predicted from quenched chiral perturbation theory. However, in the 
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'"■bare"' 










U.U l Oil 1J 1 


0.01 


0.0934(90) 


0.0889(12) 


0.02 


0.1012(53) 


0.1003(12) 


0.04 


0.1060(34) 


0.1149(12) 


0.06 


0.1106(25) 


0.1273(13) 


0.08 


0.1157(22) 


0.1389(13) 


0.10 


0.1209(21) 


0.1502(13) 



Table 6: Pseudoscalar decay constants from the "indirect" method, Eqs. M9\) . 

case of Wtm fermions we observe a bending of /* m when the pion mass is small. It 
has been argued |1] that the condition 

m q a > a 2 A 2 QCD (21) 



(where m q = yjfh 2 + \i 2 is equal to fi at full twist) has to be satisfied (with some 
proportionality factor in front of the r.h.s.) in order the explicit breaking of the 
chiral symmetry to be driven by the the mass term rather than by the Wilson term 
(in which case large cut-off effects could appear). If this is the effect we are seeing in 
Figure |Hl then this inequality seems to be satisfied with a proportionality constant of 
0(1). In order to visualize this, we represent the squared pion mass corresponding 
to a bare quark mass rh q a = a 2 Aq CD as the vertical dotted line in Figure |H1 What is 
puzzling is that, for a setup where the correlators are 0(a) improved (as it should 
be in our simulations of Wilson twisted mass fermions), one could hope that the 
condition that has to be satisfied is m q a ^> a 3 AQ CD rather than the condition in 
Eq. (|21|). In this case, all of our data should be safe, unless the proportionality 
constants in front of the r.h.s. is a large number. Clearly the bending phenomenon 
observed in Figure |H1 deserves further analytical and numerical investigations and can 
be clarified presumably only when results at smaller lattice spacings are available. 

For overlap fermions, the renormalization constant can be reliably extracted in 
the chiral limit and hence also the "direct" method can be used. For these two 
methods we have considered both the simplest correlators Cp. ov (xo) ("indirect") 
and C J 4 0>ov (xo) ("direct"), as well as the composite ones Cp-s,ov(%o) ("indirect") 
and C J 4 0+ y 0i ov(^o) ("direct"), where the finite volume effects from the zero modes 
cancel. The results in the chiral limit for the four cases turn out to be compatible 
within the errors, 

hm J™ s a = 0.0963(52) , 

m ov — »0 ° 

lim fi?a = 0.0980(46) , 

m ov ^0 

lim W a = 0.1010(40) , 
lim f^ A a = 0.1018(38) . 

m ov — >0 o 



16 







1 ' 1 ' 1 ' 1 

_CLO \ ° P V ^ ^ Q © 


— r 
\ 


c 

J 


f K a-fa(1 + o^M K 2 /(47cf) 2 ) 

1 , i , i , i 



5 I L_i I i I i I i l_ 

0.05 0.1 0.15 0.2 



(M K a) 2 



Figure 9: Difference between the results for the decay constant and a linear fit to 
degenerate data only. 

This means that the finite volume effects, which in the case of Cp }OV (xq) and Ca , ov {xo) 
may affect both the pseudoscalar mass and the matrix element, cancel out when one 
takes the suitable ratio needed to compute /„-. The slight difference between the 
"direct" and the "indirect" method can be ascribed to the different 0(a 2 ) lattice 
artifacts of the two correlators used and also to the uncertainty in the chiral extra- 
polation of Z°a needed for the "direct" method. 

Since for overlap fermions the behavior of f° v is perfectly linear, as predicted by 
quenched chiral perturbation theory — in contrast to the Wtm case — we can hope 
to check the prediction of qxPT in the case of f^. Moreover, since the two methods 
for extracting f n are in good agreement, in the non-degenerate case we use only the 
PCAC method, 

fK=^^^\(0\P\KU. (22) 
At one loop in qxPT, fx takes the form 

/* = /(l + ^M*+FV + LG) (23) 

where / and a$ are the same as in Eq. (|2Uj). "FV" are finite size effects and "LG" 
logarithmic corrections (see Ref. for the complete formula). 

One could envisage the following strategy: determine / and from the degen- 
erate data and search for FV and LG. Unfortunately, within our statistical accuracy 
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we see only a linear behavior in M^, with as in perfect agreement with the determi- 
nation from the degenerate data (i.e. from /„.). This is displayed in Figure EE where 
we plot the difference between the data and the values obtained from the q%PT 
formula for fx without the terms FV and LG and where as and / are determined 
from the degenerate case only. 

The results we get are f„ = 155(11) MeV, f K = 173(8) MeV, f K /f v = 1.11(3), 
as = 1.85(30). As usual in the quenched approximation, the value of fx/ fn turns 
out to be about 10% smaller than its experimental value. 

3.5 Baryon masses 

In order to extract baryon masses we use the following interpolating operators (for 
the octet and the decuplet respectively), 

B? = e ABC [((d A fC l5 u B )u c a - ((u A ) T C l5 d B )u°] , 

Bt = e ABC {{u A ) T C lk u B )u c a , fc = l,2,3, (24) 

where C is the charge conjugation matrix. For the decuplet k = 1, 2, 3 are equivalent. 
We have chosen k = 1. 

For correlators at zero momentum in the overlap case, we have 

+ ^U(B° a ct ' dcc (x, x Q )B°; tAcc (0)) ov oc e~ Mx ° , a « x < - . (25) 

X 

In this case we perform a simple exponential fit in the first half of the lattice in the 
time direction to avoid contaminations coming from the state with opposite parity. 
In the twisted mass case (at twist angle uj = vr/2) it is easy to show that 

(5 oct,dcc (x)i? oct,dec (0))phys = 1 (1 + Z75 ) Q7 (5f' dcC (x)i?r' dCC (0))t m (l + ijtU , (26) 

where (-B° ct ' dec (x)S^ ct ' dec (0))phys is the correlator with the correct quantum numbers 
in the continuum. 

Results (obtained by using sink-smearing) are reported in Table and plotted 
in Figure Concerning overlap fermions, due to the relatively small volume, the 
decuplet channel as well as the octet correlators corresponding to the lowest two 
bare quark masses are too noisy and we are not able to extract the corresponding 
masses. In the Wtm case, we observe a bending of the data at small quark masses 
similar to the case of the vector meson, whereas at bare quark masses larger than 
50 MeV the results seems to behave very similarly to the overlap data. 

4 Conclusions 

In this paper we confronted quenched simulation results of overlap fermions and 
Wtm fermions in their approach to the chiral limit. We emphasize that we tested 
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Figure 10: Comparison of results for the baryon masses (octet and decuplet in the 
SU(3) symmetric limit) as a function of the pion mass squared for overlap and 
twisted mass fermions. 



both lattice discretizations of QCD at only one value of (3 = 5.85 corresponding to 
a value of the lattice spacing of a = 0.123 fm. Although scaling tests are certainly 
of great importance to further explore the potential of both kind of fermions, our 
results show already very interesting features. 

The first is that indeed with both kind of lattice fermions small values of the 
light meson masses can be reached, such as M n ~ 220 MeV in the overlap case 
and M n ~ 270 MeV in the Wtm case. In addition, the statistical fluctuations for 
the observables studied here are comparable for both formulations. This is very 
promising. In a detailed algorithmic study [§] we find that Wilson twisted mass 
fermions are a factor of 20 to 40 faster than overlap fermions. Thus Wilson twisted 
mass fermions have the potential for dynamical fermions simulations at realistically 
small quark masses on the next generation of supercomputers in the multi-teraflops 
range. 

However, we believe that a number of questions have to be addressed to under- 
stand better the Wtm formulation of lattice QCD: In the quenched comparison of 
overlap and Wilson twisted mass fermions, we encountered a "bending" phenomenon 
for Wtm fermions. This effect manifests itself in the chiral approach of all the quan- 
tities studied here. While the data for overlap fermions extrapolate nicely linearly, 
the data for Wtm fermions show a bending when the quark mass assumes too small 
values. This effect might be explained by the interplay of the Wilson term and the 
twisted mass term, which requires one of the two inequalities discussed in Sec. 13.41 
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M. 



dccuplct 



a 



0.01 
0.02 
0.04 
0.06 
0.08 
0.10 



0.8479(58) 
0.8885(31) 
0.9336(19) 
0.9788(14) 



0.592(37) 

0.782(25) 

0.963(11) 

1.0987(88) 

1.2183(74) 

1.3302(72) 



0.941(65) 
1.029(36) 
1.126(25) 
1.2199(15) 
1.3124(12) 
1.4136(83) 



Table 7: Baryon masses with sink smearing. 



to be satisfied. 

For the present simulations at only one value of the lattice spacing we are not 
able to determine the cause of this "bending" phenomenon and whether it disappears 
in the continuum limit. For this, a detailed scaling analysis would be necessary, a 
work that is in progress. Overlap fermions on the other hand nicely approach the 
chiral limit close to the physical point with realistic light meson masses. It seems 
that the conceptual virtues of this approach become more and more important, the 
smaller the quark mass is chosen. However, we believe that a final answer, which 
formulation to use for extracting physics in the chiral limit, can only be given when 
the scaling behavior of both approaches is understood. 

For dynamical simulations the presence of a first order phase transition has been 
seen for Wtm fermions ^Uj- This observation is in accordance with an effective 
potential picture [211 12E1 12H1 1221 • The findings in Ref. [TU] just emphasize the fact 
that the study of the phase structure of lattice QCD is a necessary prerequisite for 
reliable physics results. 
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